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Abstract. - We identify an optical signature for detecting entanglement in experimental 
nanostructure systems comprising coupled excitonic qubits. This signature owes its strength to 
non-Markovtan dynamical effects in the second-order temporal coherence function of the emit- 
ted radiation. We calculate autocorrelation and cross-correlation functions for both selective 
and collective light excitation, and prove that the coherence properties of the emitted light do 
indeed carry information about the entanglement of the initial multi-qubit state. We also show 
that this signature can survive in the presence of a noisy environment. 

Essentially all the proposed schemes for quantum information processing - including quan- 
tum computing and teleportation - require the reliable generation of entangled states between 
pairs of qubits. Many physical implementations of qubits have been proposed, and even built, 
in systems ranging from artificial nanostructures through to naturally-occurring molecules 
and even biological systems. However the following common problem faces all such systems: 
Having gone through the effort of forming a set of N > 2 qubits with the intention of un- 
dertaking some form of quantum information processing, how can we be sure that entangled 
states are indeed being generated experimentally! More specifically, are there any signatures 
for two-qubit systems that can distinguish between states such as mixed or product states 
(which exhibit at most classical correlations) and the crucially important entangled states 
which carry purely quantum correlations? This question provides the first motivation for the 
present work. The second motivation stems from the fact that the preparation and detection 
of highly-correlated quantum states are difficult to perform in a controlled way due to the in- 
teraction of the quantum system with the noisy environment - after all, this is still one of the 
main practical hurdles facing quantum computation. There is, therefore, a general need for 
detailed quantitative theory concerning the time-evolution of different initial quantum states, 
including their decoherence properties. 

In the present work, we propose such a signature, whose experimental implementation 
can be achieved using current ultrafast optical spectroscopy. In addition to this important 
practical finding, our work provides a fundamental example of a quantum system where the 
Markov approximation, which is invariably employed in such calculations, yields incorrect 
results. Indeed, our results suggest that any such optical probe of entangled states can only 
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be correctly interpreted using non-Markovian theories. The Markov approximation, which is 
so convenient operationally, is simply not good enough in these systems. 

For the sake of concreteness, the qubit-matter systems we are addressing can be realized 
using localized excitons in semiconductor quantum dots (QDs) [1-3] or some organic/biological 
molecules [4] for which ultrafast optical spectroscopic techniques are available. Moreover, 
charge or phase superconducting qubits [5] could also be studied within the present proposal 
if the driving and detected radiation fields are in the microwave sector of the electromagnetic 
spectrum. Non-Markovian effects have been demonstrated in CuCl [6] and InGaAs[7] QD's in 
a picosecond time scale. At this short times the relaxation process cannot be described by a 
single rate constant. For example, in CuCl QDs the relaxation time is w 2ps, thus the expected 
value of the reservoir correlation time is or the order of at least lps. Additionally, ultrafast 
time-scales should be understood as femto- or picoseconds for excitons, and microseconds for 
superconducting qubits. 

The photostatistics of photon counting as measured with the correlation function g^> (T, T+ 
t) has triggered intense theoretical and experimental activity in QDs [8-11] and biological 
systems (chromophores) [12-14]. This temporal second-order correlation function gives in- 
formation about the conditional probability of detecting one photon at time T + t provided 
a previous photon was found at time T. We are especially interested in the short time dy- 
namics of such correlations where the coupled system-environment could show some unusual 
behaviour such as recoherence, among others. Therefore, a proper description of the dynamics 
must be undertaken including non-Markovian effects. Several theoretical approaches can be 
found in the literature: the unraveling of non-Markovian stochastic Schrodinger equation [15], 
non-Markovian Monte Carlo methods [16] and the time-convolutionless projection operator 
technique [17,18]. This last technique, based on the Zwanzig projection approach [19], has 
the advantage of allowing the solution of a local-in-time master equation while taking into 
account memory effects due to the system-environment coupling. We adapt this last approach 
for the coupled qubit-noisy environment which we address. We report three basic results: (i) 
Unambiguous signatures of initial multi-qubit states; (ii) detailed numerical calculations of the 
temporal auto-correlation and cross-correlation, from which new experimental features can be 
detected and analyzed at ultrafast timescales, and (Hi) the identification of induced effects 
by driving selectively an individual qubit, and its effective non-classical radiation action on a 
second nearby qubit. 

We consider the concrete example of two optically-driven, dipole-dipole interacting qubits 
in contact with a boson bath. The reduced density-operator in Lindblad form, is: 

r,. 2 2 2 2 

i=l i=l i=£j »,i=l 

where S\ is a qubit raising operator, is the Rabi frequency and A, = w^-Wi the laser detun- 
ing for the i th -qubit. -fij and Vij (i ^ j) represent the collective decay and the dipole-dipole 
interaction [20] between qubits, respectively. Using the Zwanzig [19] projection technique, it 
is possible to obtain a master equation governing the time evolution of the density operator of 
the system in which the earliers time dependence in the density operator are included in a time 
dependent decay rates. The non-Markovian effects are included by taking a time-dependent 
spontaneous decay 7i,j(t) = 7i(t) [17,19]. In Eq.JT]), where the first three terms are associated 
with the coherent evolution of the system and the last term describes decoherence processes, 
there are three control parameters: the laser field strength Pi, the driving laser detuning with 
the qubits, and the inter-qubit separation r%2 which affects Vij. We assume that the prop- 
agation direction of the incident field is perpendicular to the interqubit axis, k _L f\2. Note 
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that although our theory employs a generalization of Markovian decay-rates to non-Markovian 
situations, i.e. jij(t) = 7,(i), the strong resulting enhancement of photon-photon correlation 
effects which we observe is a highly non-trivial consequence of this generalized decay rate. 
In addition, a wide range of monotonic time-dependent functions ji(t) should show similar 
effects (but different from the 7 constant case), making our results insensitive to the precise 
form assumed for Jij(t) [21,22]. 

An important aspect of the present work is the analysis, at very short times, of the finite 
memory response of the reservoir (7^) due to the qubit-boson environment interaction [21]. 
Extending the projection technique result [8] to the present coupled qubit case, the qubit 
spontaneous decay is taken as 7j(i) = r (l — e~ lrt ), where T is the usual (stationary or 
Markov) decoherence decay constant. In order to preserve the positivity of the reduced density 
matrix, we consider an intermediate regime of coupling strengths between the environment 
and the qubit s. 

We focus on three central optical measurements: (i) the autocorrelation function (i = j), 
(ii) the cross-correlation functions (i ^ j) given by 



Physically, these functions represent the conditional probability of detecting one photon emit- 
ted from the j th qubit at time T + r after detecting one photon emitted from the i th qubit at 
T; their importance relies on the fact that they are independent of the propagation observa- 
tion angle and also that they are easier to obtain experimentally. Finally, (Hi) the full photon 
correlation is given by 



to describe the time-evolution of qubit correlations such as (Sj(T)Sj(T + T)S^(T + T)S^(T)). 
Since we describe the interacting qubit-environment dynamics as a non-Markovian process, 
it turns out that the set of equations to be solved closes in a finite hierarchy of operator 
correlations. The above correlation couples, in a nontrivial manner, with higher-order corre- 



lations of type (S^(T)Sf(T + T)S^(T + T)S^(T + T)S^(T + T)St{T)), due to the fact that 



the quantum regression theorem cannot be safely used in the present context where memory 
effects are included. 

We consider two important experimental setups: (i) a collective excitation where both 
qubits are simultaneously illuminated with the same non-resonant laser intensity, and (ii) 
a selective excitation scheme where only one qubit is driven by the laser producing non- 
classical light, which in turn can excite the second qubit or can be registered by the de- 
tector. In the first case, it is impossible to distinguish the light emitted from each qubit 
while in the second case anti-correlation effects can be observed indicating the origin of the 
detected emitted photon. Additionally, different initial qubit states are analyzed: (i) A sep- 
arate (non-entangled) state formed by the product of identical superposition states for each 
qubit, I ^qq) = |(| 0i)+ I li)) <8> (| 02}+ I I2)), where the subscripts indicate the i = 1,2 
qubit. (ii) A maximally entangled Bell state | = -4=(| 0i0 2 }+ | li^))- (hi) An cntan- 



g%(T,T + T) = 



((S](T)sI(T + t)S-(T + t)S-(T))) 
((S}(T)S-(T))(S}(T + t)S-(T + t))) 



(2) 
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Fig. 1 - The full correlation function g 2 (T,r) for a separable initial state (|^qq)) with f3 = 1.25Fo. 
(a) Non-Markovian case with j r — 2Fo, it means that the reservoir correlation time r r = h/2To = 16ps 
and (b) Markovian case (fcri2 = n). 

gled triplet state | = ^(| 0il 2 )+ | li0 2 }). (iv) An entangled singlet state (EPR-state) 
I *s> = ^(| OiU)- | li0 2 )). 

Collective excitation (f3\ = P2 — In Fig.l, the full correlation function g 2 (T,r) is 
plotted for an initial separate state I'Fqq). Both qubits are simultaneously illuminated with 
a laser excitation of 1.25T0 (Fo=20 fieV). It is possible to estimate the dipole moment by 
fiEiTTp/h = A, where Eq is the radiation electric field, t p w 6ps is the window time of the 
incident pulse and A is the pulse area of 0.247T [23], which corresponds to a typical experimental 
excitation (with /i«8 Debyes) [2]. The central laser frequency is detuned Ai = — 0.5Fo = — A 
and A 2 = 0.5Fo = A, from each qubit. A strong antibunching of the emitted radiation is 
evident at very short times. Note that, for < T < Iq , the second-order correlation function 
is rather immune to decoherence effects in the non-Markovian case (Fig. [TJ(a)) as compared 
with the Markovian case (Fig. [TJ (b)) . This is a consequence of the recoherence processes 
between the environment and the two qubits. The system decoheres, interacting with the 
bath, and at tTq « 1 a correlation revival arises. This strong antibunching effect is a non- 
Markovian process. Thus, the quantum correlations at time T + r (a hrst photon detected 
at T) are affected by the type of initial state at time T — 0. This effect is enhanced for 
individual qubit superposition states, and is a consequence of the vanishing values taken by 
the intensity correlations J^i miSjiT + ^S^T + r)). We emphasize that this effect is stronger 
for two interacting qubits as a function of r, as compared with the single qubit case previously 
reported in Ref. [8] and persists for longer T times. This can be understood by analyzing the 
initial state as: | ^qq) — ^ ( 1 1 1 , F 2 ) + |0i , O2) + |0i, F 2 ) + |li, O2))- This state comprises 
the superposition of one exciton in each qubit, one non-radiative state (zero excitons) and 
the triplet state (which is associated to a superradiant state). The contribution from the 
triplet state provides a distinctive signature which is significantly enhanced by non-Markovian 
effects. For long times (TTq = 2) an additional peak appears in the second-order correlation 
function, as a consequence of the typical oscillations given by the external parameters at 
Tr = 4(A 2 + (2/3) 2 + (Vi 2 - To) 2 )- 1 / 2 . 

In order to stress the quite different optical signatures for different possible initial states, 
Fig-EIa) focuses on an initial Bell state: | <Fs). This state does not contain the contribution 
of the triplet state which explains why the intensity of the corresponding correlation func- 
tion is about three times smaller than the one corresponding to the separate state | ^qq)- 
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Fig. 2 - The full correlation function g 2 (T,r) for an initial Bell state, (a) Non-Markovian case and 
(b) Markovian case. Parameters as in Fig. Q but now with a common Rabi frequency f3 = 2.25Fo. 

However, the non-Markovian effects can still be observed, when compared with the results 
of Fig. EJJb) for the Markovian case. The origin of this strong antibunching signature can be 
understood by inspecting Fig. [31 with the same parameters as Fig. [1] and [2J which shows the 
results for the triplet state: | ^t}- The g 2 (T, r) results yield a huge (one order of magnitude 
higher) antibunching effect at very short T times in clear contrast to the Bell state. The 
oscillations remain in the non-Markovian case, as compared with the Markovian case (Fig. 
EJb)). Our results clearly show one distinctive signature: the triplet state is responsible for 
a strong antibunching effect, which is enhanced due to the recoherence processes between 
the environment and the qubit system. This effect can be tailored and also enhanced if the 
inter-qubit separation is reduced. We emphasize that our parameter values are consistent 
with current experimental data, and the results presented can be produced with low laser 
excitation provided that (3 > j r . 

Summarizing, at very short times the slope of g 2 (T, r) is positive indicating an antibunch- 
ing behaviour, which is enhanced by non-Markovian effects in the range of < tTq < 1 and 
< TTq < 1, allowing the distinguishability among the | ^qq},\ ^b) and | ^>t) states. As 
can be seen in Fig [3] the triplet state lacks robustness to decoherencc effects. The EPR-statc 
(| ^s)), does not show significative differences in its time-evolution when memory effects are 
taking into account (not shown) . It is also interesting that the second-order time correlations 
reach higher values for separable states but decrease when the input state becomes entangled. 
Finally, for TTo > 2 the non-Markovian results tend toward the Markovian ones, emphasizing 
that in this long time regime the time-correlation is stationary and independent of the initial 
states. 

Selective excitation (f3\ = (3,02 = 0).- We now consider the case when only one qubit 
is externally driven by a laser light and each qubit is out of resonance. In this scheme, it 
is possible to obtain information about the signatures of the light emitted from a second 
qubit when antibunched light emitted from the excited qubit is driving the second qubit. 
In Fig. ID we show comparative results for different initial states (TTq — 1), where the non- 
Markovian signatures can still be observed. In Fig.|4ja) and (b), the autocorrelation functions 
are completely different, due to the fact that the non-illuminated qubit does not emit photons. 

(2) (2) 

Once the non-driven qubit becomes excited, g\i (r) decreases while g\^ (r) increases. However 

it is interesting to note that both autocorrelations vanish at r = 0, but start to increase towards 

(2) (2) 

unity for long times. This leads us to consider the cross-correlations g\ 2 ( T ) an d 521 ( T ) j which 
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allow to know if the information about the emitted correlated photons and if the emitted 

(2) 

photons are uncorrelated {g\-{r) < 1) or correlated (for i ^ j). The differences are clearly 




2.5 




Fig. 3 - The complete correlation function g 2 (T, r) in the vertical axis for the excitonic triplet state, 
(a) Non-Markovian case and (b) Markovian case. Antibunching is comparable with a superposition 
state. 



visible in Fig. 4(c). For separate state (| 'Fqq}) the plot shows not difference between gi2 
and g2i which correspond to a complete uncorrelated situation for any r. The | ^b) state, 
exhibits strong correlations for tTq < 0.5 (312 > 2) with a time dependent correlations going 
to zero. However, by selectively exciting just one qubit, this strong correlation remains for 
longer times (up to tTq ps 2, see Fig. 4(d)). The anticorrelation is perfect for the entangled 
triplet and singlet states, indicating that there are no photons available to excite the second 
qubit and thus it is possible to distinguish the source of the detected photons. This effect can 
also be observed in the stationary regime (T — > 00), as is clearly seen in the inset of Fig. 4(a) 
for times tTq < 0.5. However, the memory is lost due to the interaction with the environment 
and therefore no information about the initial state is preserved. The cross-correlation, 921 
goes to zero at a very short time (tTq 0.2.) in contrast with the non-Markovain case where 
for a long period of time, the cross-correlations remain almost at zero value for the entangled 
singlet and triplet qubit states. 

To summarize, we have shown that for an optically-driven coupled qubit system, differently 
prepared initial states - for example, separate superpositions, Bell and triplet (symmetric) 
states - can be clearly distinguished as a direct result of non-Markovian effects. The cal- 
culation of normalized autocorrelation and cross-correlation second-order functions, also give 
valuable information concerning the extent to which each qubit behaves like a single quantum 
light source. We have demonstrated how it is possible to assess the entanglement content by 
observing the correlations function of the emitted light for collective as well as selective exci- 
tations. In the latter case, the cross correlations give enough information about the photon 
source by exploiting the non-classical light produced by a qubit on a nearby second single 
qubit. This highly correlated light could be used to probe systems sensitive to non-Markovian 
effects. Our study exposes new phenomena that can be observable at ultrafast scales of time 
for which non-Markovian dynamics should be dominant. To our knowledge, the proposed 
signature is the first of its kind and should act as a powerful tool for experimental groups 
attempting to realize quantum information processing in a wide range of optically-driven, 
coupled-qubit systems. 

* * * 
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Fig. 4 - (on-line color). Selective escitation with a Rabi frequency of fj\ = 2.25Fq. Interqubit 
separation as in Fig. ([TJ. (a) and (b) show the auto-correlated functions and (c) and (d) are the 
cross-correlated ones. The different signals correspond to a superposition (black-line), Bell (red-line), 
symmetric (green-line), anti-symmetric (blue) states. The inset depicts the stationary limit. 
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